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ABSTRACT 

Since the suggestion of relativistic shocks as the origin of gamma-ray bursts (GRBs) 
in early 90's, the mathematical formulation of this process has stayed at phenomeno- 
logical level. One of the reasons for the slow development of theoretical works has 
been the simple power-law behaviour of the afterglows hours or days after the prompt 
gamma-ray emission. It was believed that they could be explained with these for- 
mulations. Nowadays with the launch of the Swift satellite and implementation of 
robotic ground follow-ups, gamma-ray bursts and their afterglow can be observed in 
multi-wavelength from a few tens of seconds after trigger onward. These observations 
have leaded to the discovery of features unexplainable by the simple formulation of 
the shocks and emission processes used up to now. Some of these features can be 
inherent to the nature and activities of the GRBs central engines which are not yet 
well understood. On the other hand devil is in details and others may be explained 
with a more detailed formulation of these phenomena and without adhoc addition of 
new processes. Such a formulation is the goal of this work. We present a consistent 
formulation of the kinematic and dynamics of the collision between two spherical rel- 
ativistic shells, their energy dissipation, and their coalescence. It can be applied to 
both internal and external shocks. Notably, we propose two phenomenological models 
for the evolution of the emitting region during the collision. One of these models is 
more suitable for the prompt /internal shocks and late external shocks, and the other 
for the afterglow/external collisions as well as the onset of internal shocks. We cal- 
culate a number of observables such as flux, lag between energy bands, and hardness 
ratios. One of our aims has been a formulation enough complex to include the essential 
processes, but enough simple such that the data can be directly compared with the 
theory to extract the value and evolution of physical quantities. To accomplish this 
goal, we also suggest a procedure for extracting parameters of the model from data. 
In a following paper we numerically calculate the evolution of some simulated models 
and compare their features with the properties of the observed gamma-ray bursts. 
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1 INTRODUCTION 

The Swift (|Gehrels et all [2004) observations of more 
than 200 Gamma-Ray Bursts (GRBs) and their follow- 
ups have been a revolution in our knowledge and under- 
standing of these elusive phenomena. The rapid slew of 
the S wift X-ray and UV/ optical teles copes - respectively 
XRT (|Burrows et ql.ll2005h and UVOT (|Roming et 01.1120051 ) 
- as well as ground based robotic telescopes have permitted 
to observe GRBs and their afterglow in multi-wavelength 
from few tens of seconds after the pro mpt or precursor 
gamma-ray emission is detected by BAT (|Barthelmv et all 
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2005), up to days after trigger. These observations show that 
the emission can be essentially divided to three regimes: 1) 
The prompt gamma-ray emission which can be very short, 
few tenth of milliseconds, or long, up to few hundred of 
seconds; 2) A tail emission in X-ray which is observed for 
more than 90% of bursts. For some bursts this tail is also 
detected as a soft faint continuum in gamma-ray. In about 
40% of bursts this early emission has been detected in op- 
tical and infrared too. In this regime for many bursts flares 
have been observed mainly in X-ray. Sometimes the counter- 
part of flares have been also observed in gamma-ray and/or 
optical/IR. In many bursts the early steep slope of the X-ray 
emission at the beginning of this regime becomes much shal- 
lower and somehow harder at the end; 3) The late emission 
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can be considered as the epoch after the break of shallow 
regime in which the emission is usually a continuum and 
no or little flaring activity is observed (but there are excep- 
tions su ch as GRB 070110 |Krimm et all 120071 ') and GRB 
081028 (|Guidorzi et ffll.ll2008t ) which had bright late flares). 
The duration and relative fluxes of these regimes can vary 
significantly between GRBs. 

In one hand, it seems that the idea of synchrotron 
emission from accelerated particles in a relativistic shock 
as the orig in of the prompt emission l|Rees fc Meszarosl 
1 1994 [ 1998) is essentially correct. On the other hand, 
the early observations of what is usually called the af- 
terglow - the emission in lower energy bands usually ob- 
served from < 100 sec after trigger onward - have been 
the source of surprises and raised a number of que stions 
about many issues: the activity llFan fc Weil l2005h and 
the nature of the engine (|Petrovic et all 120051 ; iFrver et all 
l200d : iGal-Yam et J z1 l2008l ) , the concept of prompt /interna l- 
afterglow/external shocks l| Ramirez- Ruiz fc Granotl l2007h . 
the efficiency of energy tran sfer from the outflo w -the fireball 
- to synchrotron ra diation |Zhang et aZ.|[2006h . the collima- 
tion and jet break (|Covino et 01.1120061) . the behaviour of X - 
ray and optical light curves ( Panaitescu fc Vestrand| [2008). 
etc. Many of predictions such as the existence of a signifi- 
cant high latitude emission with a strict relation between the 
light curve time evolution slope and the spectrum index, and 
an achromatic jet break have not been observed. Moreover, 
the origin of unexpected behaviours such as a very steep 
decline in low energy bands after the prompt (|Zhang et al\ 
l2007bl ) and a very shallow regime which lasts for thousands 
of seconds are not well understood. Other unexpected ob- 
servations are the existence of the chromatic multiple breaks 
in the X-ray light curves, flares in X-ray and optical bands 
hundred of seconds after the prompt even in some short 
burst s (ex: GRB 0603 13 llPagani et all [20061: [ Roming et all 
l2006h . GRB 070724A (|Ziaeepour et all 2007all. a tail emisT 
sion in short b ursts (ex: GR B070714B (IRacusin et aZ.ll2007l ; 
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20081), GRB0804 26 (|Ziaeepour et all l2008bh . GRB 



080503 (|Mao et aZj |2008)). and very short, hard, and high 
amplitude spikes in long bursts that could lead to the clas- 
sification of the burst as short if the instrument was not 
enough sensitive to d etect the rest of the prompt emis- 
sion (ex: GRB06061 4 dMangano et qZj|2007l : iGehrels et al\ 
l2006h . GRB061006 (|Schadv et aZj|2006ah ). This makes the 
classi fication of bursts as short and long much more ambigu- 
ous (|Zhang et oZ.||2007af ). 

One conclusion that has been made from these obser- 
vations is that the central engine can be activ e for up to 
thou sands of seconds after the prompt emission (|Fan fc Weil 
2005). But the nature of the fireball and its source of en- 
ergy is not yet well understood, and we can not yet ver- 
ify this interpretation or relate it to any specific process 
in the engine. It seems however that whatever the ori- 
gin of the fireball, it must be baryon dominated other- 
wise it could not make long term effects correlated to the 
prompt emission. In this case, the internal and external 
shock models as the origin of the prompt and afterglow 
are good candidates. Nonetheless, the lack of a simple ex- 
planation for the observed complexities has encouraged au- 
thors to consider other possibilities, for instance associat- 
ing both the prompt gamma-ray and the afterglows to ex- 
ternal shocks and fast variations to abrupt density varia- 



tion of the surrounding material (Ra mirez-Ruiz fc Granotl 
l2007fl . However, it has been shown that in such models it 
is not possible to explain the fast variations of the prompt 
even in pres ence of a bubbly environment or pulse-like den- 
sity c hange dRamirez-Ruiz fc Granotll2007l : lNakar fc Granotl 
12007ft . 

Here we suggest that at least some of the features 
of early afterglows can be related to a complex shock 
physics and/or features in the fireball/jet. In fact, sim- 
ulations of the acceleration of electrons and positrons 
by the first and second Fermi processes show that the 
evolution of electric and magnetic fields as well as 
the energ y distribution of acce l erated particles are quite 
complex dBednarz fc Ostrowskil Il996l : iRieger et all 120071 : 
iDieckmann et aZ.ll2006h . Plasma instabilities lead to the for- 
mation of cohere nt electric and magnetic fields and acceler- 
ation of particles ( Yang et aZ.lll994 | ; | Wiersma fc Achterberel ; 
ISilva et aZll2002l : lReville et aZ.ll2006h Their time evolution 
in relativistic shocks can significantly affect the behaviour 
of the prompt and the afterglow of GRBs. If the number 
density of particles in the ejecta is significant and the shock 
is collisional, the state of matter in the jet can be also an 
important factor in determining the behaviour of the fields, 
and thereby the synchrotron emission by accelerated elec- 
trons and positrons. Many aspects of these processes are 
not well understood, however realistic interpretations of ob- 
servations should consider these complexities at least phe- 
nomenologically. For instance, the simple distributions such 
as a power-law distribution for Lorentz factor of electrons, or 
a constant magnetic field for the whole duration of prompt 
and afterglow can be quite unrealistic. Ideally, these quanti- 
ties should come from the simulation of Fermi pr ocesses and 
plasma instabilities such as W eibel instability (|Yang et all 
1994; W iersma fc Achterberel) that produces the coherent 
transverse magnetic field. However, these phenomena are 
complex and their simulations are very time and CPU con- 
suming. For these reasons they can not yet explore the pa- 
rameter space of the phenomena and are mostly useful for 
demonstrating the concepts and how they work. Therefore 
we are obliged to use simple analytical approximations for 
quantities related to the physics of relativistic shocks. In this 
situation a compromise between complex non-analytical ex- 
pressions and too simplistic and too simplified but unreal- 
istic analytical behaviour of the physical quantities can be 
the consideration of intervals in which a simple analytical 
function can be a good approximation. Then, by adding to- 
gether these intervals - regimes - one can reconstruct the 
entire evolution of a burst and its afterglow. 

Even with a simplified presentation of the physical 
processes one would not be able to explain GRB data 
without a model including both microphysics and dynam- 
ics of the fireball. The majority of works on the mod- 
elling of shocks an d synchrotron emi ss ions eithe r dea l 
with the emission JSari fc Piranl 1 19951: ISari et all Il996l . 



1 19981 : iNakar fc Piranl 120041 : IZhang fc Kobavashil 120051) or 

with the kinemat ics of the shock dBlandford fc McKeel 

1 19771 : iPiran et all 1 19931 : iFenimore et all 1 19961), o r both 
but i n a phenomenological way dRhoada 1 19990 ■ Few 
works (|Daigne fc Mochkovitch1ll998l . 120031 ') have tried to in- 
clude both these aspects in a consistent model, but either 
they have not been very successful - their predictions spe- 
cially for quantities such as lags in different bands were far 
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from observed values and additional parametrization was 
necessary - or the formu lation is too abstract to be com - 
pared directly with data l|Vietrill2003l : iBlasi fc Vietrill2005h . 

With these issues in mind, in Sec [5] we present a sim- 
plified shock model that includes both the kinematics of 
the ejecta and the dynamics of the synchrotron emission. 
The microphysics is included by the means of a simple 
parametrization. We calculate a number of observables such 
as flux, hardness ratios, and lags between different energy 
bands. In this paper and paper II in which we simulate part 
of prompt and afterglow regimes of GRBs in some time in- 
tervals, we show that this model can explain many aspects of 
bursts as long as we divide the data to separate regimes. The 
reason is that the simple parametrization of microphysics 
in this model can be valid at most in a limited time in- 
terval. Each regime should be separately compared to ana- 
lytical and numerical results for extracting the parameters. 
The results will show how parameters that are considered 
as constant in this model evolve during the lifetime of the 
burst. This is the best we can do until a better understand- 
ing of relativistic shock models and Fermi processes become 
available. If the model and the estimation of its parameters 
for each regime is sufficiently correct, adding them together 
should give us an overall consistent picture of characteris- 
tics of the burst, its afterglow, and its surrounding material. 
Apriori this knowledge should help to better understand the 
engine activities and eventually its nature and classification. 

The model presented here depends on a large number 
of parameters and we need an extraction procedure permit- 
ting to extract as much as possible information about the 
physical properties of the shock from the available data. In 
Sec[3]we explain how in the frame work of this model one 
can extract various quantities from data. Evidently the suc- 
cess of the modelling strongly depends on availability and 
quality of simultaneous multi-wavelength observations. 

The Swift observations show that during the first few 
hundred seconds after the trigger there is usually a very 
close relation between the prompt gam ma-ray emission 
and the emission in lower energy bands (Meszaros 2007; 
iKumar et all 120061 ) . therefore, most probably they have a 
common origin, presumably internal shocks. However, his- 
torically and even in the present literature (and sometimes 
believes) any emission after the prompt gamma-ray is called 
the afterglow - meaning due to a shock with ISM or sur- 
rounding material, presumably external shocks. Therefore, 
for clarity of context here we define the afterglow as the 
emission in any energy band and at any time after the main 
prompt peaks regardless of its origin. If by afterglow we mean 
the external shocks, this is mentioned explicitly in the text. 

We finish this paper with some outlines and two appen- 
dices containing the details of calculation of the dynamics 
and flux for power-law distribution of electrons Lorentz fac- 
tor. 



2 SHOCK MODEL 

In this section we first give a sketch description of a rela- 
tivistic collision between two shells of material and processes 
which produce gamma-ray and other radiations. Then, we 
discuss a simplified mathematical formulation of the evolu- 
tion of the shock and synchrotron emission. By restricting 



the model to a thin layer and to the early times after begin- 
ning of the collision, we can analytically calculate various 
observables. 



2.1 Qualitative description of a relativistic shock 
and its simplified model 

We begin this section by a pictorial description of present 
believes about the origin of GRBs and their afterglow. A 
central source - supernova, collapsar, collision of two com- 
pact objects, etc. - ejects highly relativistic cold baryon 
dominated material as a spherical, jet or torus-like shells 
call the fireball. Collision between faster later ejected shells 
with slower earlie r ejected ones produces what is called the 
prompt emission (|Rees fc Meszajodll994h . 

Apriori, there is no reason why faster shells should be 
ejected later. One way of expla ining this paradigm is the 
decel eration of the front shells |Fenimore fc Ramirez-Ruizl 
1999) by surrounding material which are observed 
aroun d massive objects such as Wolf-Rayet (WR) 
stars (Eld rigde et all 20061 ) . It is possible that this initial 
deceleration is the source of a weak emission which has 
been seen before the main spike in many bursts. A rel- 
atively weak and soft precursor spike has been also ob- 
served hi_s£rne_l2ursts_and can be related to this deceler- 
ation (|Umeda et all 20051). Although o t her o rigins such as 



jet-star interacti on (Lazzati et al. 20061. [20071 ) and fallback 
to the collapsar (|Wang fc Meszarosll2007l ) are also suggested 
to explain precursors, deceleration of the initial shell seems 
to be a more natural explanation and does not need any fine 
tuning of the progenitor models and their parameters. In 
contrast, jet-star interaction scenario can not explain large 
time lag - few 100 of seconds between precurs or and the 
main spike in some bur sts e.g. GRB 50820A dPage et all 
| 2005l;ICenko et aZ.ll2006l ). GRB 060124 dHolland et aZ.ll2006l; 
Romano et al.ll200Sl ). and GRB 070721B (|Ziaeepour etall 
l2007bl ). In the fallback scenario a weak jet is produced dur- 
ing the formation of a transient neutron star which later 
collapses to a black hole. The main jet in this scenario is 
produced by the accretion of the material from a disk to the 
black hole. In this case, the lag depends on the lifetime of 
the proto-neutron star and the rate of the accretion from the 
surrounding disk. Although these parameters can be tuned 
to explain the lag, apriori much longer lags should be also 
possible but never observed. In the deceleration scenario the 
maximum possible lag is the duration of the central engine 
main activity - few hundred of seconds according to the ob- 
servations of the main flares in X-ray, and is consistent with 
all the observations. The UV emission from the precursor 
should ionize the unshocked material in front of the first 
shell and therefore the re would be little ad ditional absorp- 
tion of s oft X-ray later (I Watson et aZ.ll2007l). In f act in GRB 
060124 (|Holland et aZfeOOd ; iRomano et aZ.ll20Qg| ) which had 
long lag between the precursor and the main peak, a slight 
increase in Nh column density at late time with respect to 
the initial density has been detected. In Paper II we argue 
that another possible origin of the precursor is the tempo- 
rary dynamical reduction of the emission in the early stage 
of the shock that make the burst unobservable for a short 
time. Then, with the progress of coalescence of the shells the 
emission resumes and is observed as the main peak. 

When two high density shells collide, in the most gen- 
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eral case they partially coalesce. Then, the most energetic 
particles get ahead of the rest in the downstream, and at 
the end of the collision the configuration includes again two 
shells. The back shell consists of slower particles that have 
lost their kinetic energy (shocked particles). The front shell 
is the remnant of faster unscattered particles in the shells. In 
practice we expect that the kinetic energy difference be con- 
tinuous and slower particles become an expanding tail be- 
hind a faster and most probably denser head. In few b ursts 
such as GRB 060607A Jziaeepour et all l2006bl . l2008al ) and 
GRB 070107 |Stamatikos et all 120071 ) it seems that we are 
seeing the separation between these components in the X- 
ray afterglow. Fig[T] shows a sketch of the shock processes. 

In the simplest case the shock between two shells is 
radiative. This means that for an observer in the rest frame 
of the fast shell the kinetic energy of the falling particles 
from the other shell is immediately radiated and particles 
come to rest and join the shell. For a far observer at rest 
with respect to the engine the difference between the kinetic 
energy of the two shells is partially radiated and partially 
transferred to the particles of the slower shell. The fast shell 
is decelerated until the totality of the slow shell is swept. 

Not all the particles in the fast shell are decelerated with 
the same rate. Therefore after the coalescence of the shells 
there is a gradient of Lorentz factor from the shell front 
head with highest F to the back tail with lowest F. In the 
macroscopic treatment of the shock processes it is usually as- 
sumed that during the collision two distinguishable shocked 
layers and corresponding discontinuities are formed in each 
side of the boundary between the shells. They are called 
forward and reverse shocks according to their ev olution di- 
rection with respect to the initial discontinuity dSari et all 
1 19961 ; iNakar fc Piranll2004 IZhang fc Kobavashil I2005T I . Ac- 
cording to these models and depending on the density dif- 
ference of the shells and their relative Lorentz factor, the 
induced electric and magnetic fields and thereby the distri- 
bution of accelerated electrons and their synchrotron emis- 
sion in these shocked regions can be very different. Notably 
the reveres shock is expected to emit mostly in optical wave- 
length with a relatively large lag with respect to the prompt 
gamma-ray as the emission must traverse the width of both 
shocked layers. Therefore it should app ear as an optical flash 
asynchronous from gamma-ray peaks (|Kobavashill200oh . 

On the other hand, if shells have similar densities and 
a small relative Lorentz factor, the scattering of particles 
in the two sides of the boundary quickly homogenizes the 
shocked regions, and therefore one can consider a single 
shocked zone with a relatively slow gradient in density and 
fields. If one of the shells has a density much larger than the 
other, the width of its shocked layer would be very small with 
respect to the shocked region in the other shell and again the 
assumption of a single shocked layer is a good approximation 
for internal shocks. This schematic view corresponds well to 
what is observed in the simulations of Fe rmi accelera t ion in 
the ultra-relativist i c shoc ks (see Fig.l of ISpitkovskvl (|2008T ) 
and iKeshet et all (|2008l 'l1. 

For a far observer who only detects photons from the 
synchrotron emission of the shocked material it is very dif- 
ficult to distinguish between photons coming from distinct 
regions unless they are well separated in time and in en- 
ergy band. The lack of a clear evidence for a reverse shock 
emission in the Swift bursts, specially during the prompt 



emission, means that the reverse shock in the prompt GRB 
emission is weak and the assumption of just one shocked 
region is a good approximation. Most of the synchrotron 
emission is expected to be emitted by charged particles in 
the shocked region. But as electric and magnetic fields, and 
accelerated particles are not really confined in this region 
and penetrate to a larger area, the region that emit radi- 
ation can be much extended than shocked region. We call 
this emitting zone the active region. 

Evidently, in a real situation multiple shells are ejected 
in a short time interval. In this case both pair collisions be- 
tween late shells and collision- coalescence of the late shells 
with the main ejecta - prompt shell - is possible. These 
events happen at different points of space-time, but their 
synchrotron emission can arrive to the observer separately, 
partially overlapping, or simultaneously. Therefore it is not 
always possible to distinguish separate collisions and their 
characteristics. To this complexity one must also add the 
variation of quantities such as density and Lorentz factor 
in a single shell. They increase the variability of observed 
emissions. On the other hand, overlapping emissions make 
the comparison of data with the models based on a simple 
parametrization of the physical properties of a shell ambigu- 
ous. Despite these complexities one should be able to con- 
sider peaks as separate collisions and find an effective set of 
parameters to model each one separately. 



2.2 Shock evolution 

A shock is defined mathematically as a discontinuity in 
the density distribution of a flow. It should satisfy at each 
point of the space-time the total energy- momentum a nd cur- 
rent/particle flux conservation equation (|Anildl 19891 ): 
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where T^ v is the total energy-momentum tensor; p is density 
and velocity vector; the index i indicates species of par- 
ticles/fluids with a conserved number. They also include the 
interactions between these particles/fluids. When there is no 
interaction, these equations must be satisfied separately for 
each species. In particular, at the shock front where these 
quantities are discontinuous the conservation equations take 
the form of jump conditions across the discontinuity surface: 
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The symbol [ ] means the difference of the quantities be- 
tween square brackets on the two sides of the shock front. 
Solutions of these equations determine the evolution of the 
shock front. As for the state of the shocked material be- 
hind the discontinuity, when there is no energy dissipation 
jump conditions can be used to obtain a self-similar solu- 
tion l|Blandford fc McKed li"977T l . In presence of energy in- 
jection or dissipation however the self similarity solutions 
are only approximations and in general an exact self simi- 
lar solution does not exist. Moreover, energy dissipation by 
synchrotron emission changes the form of the conservation 
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Figure 1. Sketch of a relativistic shock: a) Beginning of the shock: A high Lorentz factor cold shell (violet) moving from left to right 
collides with a slower shell or ISM (blue), b) During the passage of the shells through each other, at the shock front an active region 
is formed where large electric and orthogonal magnetic fields are induced by plasma instabilities and charged particles are accelerated. 
They lose part of their kinetic energy as a synchrotron emission. The concept of active (emitting) region is not ne w and is commonly 
used i n the context of the formation of cosmic rays in relativistic and non-relativistic shocks by Fermi processes ifvirt anen fe Vainiol 
l2005bl lah. A fraction of these highly accelerated particles escape to the downstream. By contrast, particles that have lost their energy 
move toward upstream. This process extends the active region but reduces the gradient of quantities such as density, thus gradually 
weakens the shock, c) and d) show two possible outcome of the collision: total coalescence of the shells (radiative collision) c) and when 
after the passage of the fast shell a left-over slow shell is formed behind. In these figures colour gradient presents the velocity and density 
distributions: darker colours correspond to higher velocity/density. 



equations, i.e. it can not be treated as an inhomogeneous 
term in the differential equations describing the dynamics. 
Thus, it is not possible to extend the solutions of the non- 
dissipative case to this case even as an approximation. We 
will discuss this issue in details later in this section. There- 
fore in a dissipative shock the jump conditions are applied 
only at the initial time when a shell meets another shell or 
the ISM, and they should be considered as boundary condi- 
tions when conservation equations ((TJ) and are solved. 

Intuitively we expect that with time the discontinuous 
distributions of quantities in the shock front i.e the initial 
jump conditions change to continuous distributions which at 
far downstream distances should asymptotically approach to 
the slow shell values, and at far distances in the upstream 
direction to the characteristics of the fast shell - see Fig[JJ 
for a schematic illustration. In this transient region insta- 
bilities form electric and transversal magnetic fields. They 
accelerates electrons which subsequently lose their energy 
by synchrotron and inverse Compton emission. A far ob- 
server receives the signature of the shock mainly through the 
detection of these radiations as 7-ray burst or X-ray flare. 
Therefore finding the evolution of physics of this region is 
the main purpose for solving conservation equations. 



A full solution of equations {TJ and ([2} with the initial 
conditions ((3} and ([4]) apriori include all the necessary in- 
formation about the physical processes and their evolution. 
However, the complexity of the problem can not permit to 
solve them analytically. On the other hand, numerical simu- 
lations are both complex and time consuming, and it would 
be very difficult to cover the full parameter space and ob- 
tain results that can be compared with the observations. At 
present, simulations are only able to demonstrate the va- 
lidity of ideas about processes involved and the role of the 
instabilities in the forma tion of coherent fields (|Spitkovskvl 
I2008L iKeshet et aZ]l2008l ) . 

Here we consider an intermediate strategy. We assume 
a spherica{]] thin and optically transparent active region. 
Its average distance from the central source is r(t) and its 
thickness Ar(t). We neglect the variation of quantities in- 
side the active region and consider the average value through 
the region. This means that the evolution depends only 

1 Most of the formulation presented here is also applicable to a 
beamed ejecta/jet if the transverse dispersion of matter in the jet 
is negligible with respect to the boost in the radial direction. We 
consider the effect of beaming in Sec l2.6l 
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on r(t) the mean distance of the active region from cen- 
tral engine. As for the energy dissipation, we assume a ra- 
diative shock i.e. for an observer at rest with respect to 
the active region the incoming particles lose their energy 
through synchrotron radiation and join this region. Despite 
possibility of large contribution f rom Compton cooling of 
electrons in the GRBs promp t llStern fc Poutanenl |20o3 ; 
iPiran et al 200Sh. X-ray flar es |Kobavashi et aZj|2007h . and 
afterglow ( Wang et qZ.I l200ll ) . there is no strong evidence 
of this process in the Swift data. Cases su ch as the early 
flare of GRB 050406 dRomano et aZj|2006h that was once 
attributed to Compton cooling of a reverse shock is now 
understood to be related to the late activity of the central 
engine like all other observed flares and have the same ori- 
gin as the gamma-ray. It is also suggested that the anoma- 
lous behaviour of GRB 060218/SN 2006aj (ISoderberg et all 
2006) is due to inverse Compton cooling (|Dai et all 120061 1 ■ 
but this is an exception. Therefore here we only consider 
synchrotron radiation in a relativistic shock as the source of 
emission. 

From now on we consider the active region as an iso- 
lated shell of baryonic material. The effect of ISM/slow shell 
appears as an incoming flux, and the effect of diffusion of 
shocked material into the upstream is included in the evolu- 
tion of the thickness of the active region and other physical 
parameters such as density variation with time. The initial 
value of the density and Lorentz factor correspond to the 
values in the fast shell i.e. the shock front at the time of en- 
counter between the two shells (internal shock) or fast shell 
and the ISM (external shock) . We also assume that thermal 
energy and pressure are negligible with respect to the rel- 
ativistic boost. This assumption of cold matter is specially 
justified for the prompt emission because the temperature 
of ejected material from the progenitor is expected to be 
at most a few hundred MeV, where the kinetic energy of 
relativistic baryons must be of order of few GeV or larger. 
This approximation is not probably suitable for late inter- 
action of the shells when they have lost most of their kinetic 
energy and probably turbulence and scattering have trans- 
ferred part of this energy to thermal. 

With these approximations we write the energy- 
momentum conservation equations for the active region in 
the rest frame of the slow shell. The reason for choosing this 
frame is that the same formulation can be applied to exter- 
nal shocks where the slow shell is the ISM or surrounding 
material around the engine. The latter is considered to be 
at rest with respect to the observer (after correction for the 
expansion of the Universejj. Note that for the consistency 
and conservation of energy and momentum we have to inte- 
grate equations (fjQ) and along the active region. As we 
only consider the average value of quantities, the integration 



2 Through out this work frame-dependent quantities with a prime 
are measured with respect to the rest frame of the slow shell and 
without prime with respect to a far observer at the redshift of the 
central engine. Parameters used for parametrization do not have 
a prime even when parametrization is in the slow shell frame. 



is triviaQ: 

d(r'VAr' 7 ') ,/ l2 d(n'Ar') 

■ = 7 r — + 2r (rt Ar ) 

\ dr' 



dr 



r (n Ar ) — — = n (r)r 



dE' sy 



<hrmc 2 dr' 



(5) 



d(r' 2 n'Ar'f'f3') _ a , , , , 2 d(n'Ar') 



dr' 



a' < ( <- 



r' 2 (n'A< 



dr' 
dEL, 



+ 2r'(n'Ar')) + 



(6) 



dr' 4nmc 2 dr' 

where r' is the distance from the central engine, n' is the 
baryon number density of the fast shell measured in the 
slow shell frame, n' is the baryon number density of the 
slow shell in its rest frame and in general it can depend 
on r'. Here we assume that n' (r') = N (r' /r' )~ K . For ISM 
k = 0, i.e. no radial dependence. For a wi nd surrounding 
the c entral engine k = 2 is usually assumed l|Chevalier fc Lil 
2000). For a thin shell or jet expanding adiabatically also 
k — 2 if we neglect the transverse expansion in the case of a 
jet (collimated ejecta) . But for the collision between two thin 
shells if the duration of the collision is much smaller than 
r' /c we can neglect the density change due to expansion 
during the collision, and assume k — 0. Ar' is the thickness 
of the shocked synchrotron emitting region, 7' is the Lorentz 
factor of the fast shell with respect to the slow shell, j3' = 
\J '-y' 2 — I/7', m = m p + m e ~ m p , E' sy is the total emitted 
energy, and c is the speed of light. The evolution of the 
radius with time is: 



r'(0 - r'(t' ) = c / f3'(t")dt" 



(7) 



where the initial time t' is considered to be the beginning 
of the collision. 

Apriori we should also consider a conservation equation 
for the baryon and lepton numbers. However, with approxi- 
mations explained above the thickness of the active region is 
a parameter which is added by hand. The simplest assump- 
tion is that the active region is formed only from particles 
that fall from the slow shell to the shock front and they stay 
there. In this case the baryon number conservation is simply: 

d(n'Ar') 



dr' 



= n 



(8) 



On the other hand the assumption of particles staying for 
ever in the active region seems quite unphysical because 
scattering, acceleration, and dissipation move shocked parti- 
cles to both upstream and down stream and gradually many 
of particles are in a region where instabilities are too weak to 
make the electric and magnetic fields necessary for the accel- 
eration and synchrotron emission. Therefore in this approx- 
imation the active region can not be considered as an com- 
pletely isolated system with full conservation laws applied 
to it. The consequence of this is that we can not determine 
Ar from first principles. This point can be also interpreted 
as the manifestation of the fact that there is not an abrupt 
termination of the active zone, and therefore there is no con- 
servation for the artificial boundary we have added by hand. 

3 Note also that apriori we should consider an angular term pre- 
senting the collimation of the radial element. But for an infinites- 
imal element the angular dependence is negligible and we do not 
add it to this formulation. 
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In fact it is equally valid if we consider a conserved num- 
ber of particles and the corresponding conservation equa- 
tion but give up energy or momentum conservation. In this 
definition the active region follows active particles (in an 
Hamiltonian formulation sense). However, the only measur- 
able quantity for a far observer is the energy dissipated as 
radiation. Therefore, it is more useful to define the active re- 
gion based on the energy-momentum conservation and leave 
the number of these particles as a free parameter. Note also 
that in the left hand side of ([5]) and ([6]), as well as in the 
synchrotron term under some conditions (see below), Ar' 
always appears as n'Ar' i.e. the column density. It is more 
relevant for the observations and does not depend on the 
way we define the active region. 

We can formally integrate equation ([5jl by replacing the 
synchrotron term in the left hand side of (|6j , and determine 
the column density of the active regiorQ 



n Ar = 



1) + (3 - K)r ,2 n'(r' )Ar'(r' ) 7 ' (l - O ) 



(3 



r' 2 7 '(l - 



(9) 



where 7q = 7'(?"q) and /3' is the corresponding f3' . This solu- 
tion depends on the evolution of 7' which can be obtained by 
solving (B), but we first remind the dependence of the syn- 
chrotron term on the microphysics of the shock. The power 
of synchrotron radiation emitted by the active shell is: 



P' 



dE' sy 



c/3 



dr' 



16tt 



/2 a ' ' 

r Ar UTcy 



(10) 



where n' e is the number density of accelerated charged lep- 
tons - electrons and possibly positrons - with a Lorentz fac- 
tor 7 e 2> 7' in the active region (shock) frame, B' is the 
magnetic field - B' 2 /8tt is the magnetic energy density in 
the active region frame, and ot is Thompson cross-section. 
We define the normalization of electron distribution as the 
following: 



ra' e (7e)d7 e = 
len' e (j e )d'y e = 



'y' 2 m p n' e e 



(11) 
(12) 



where n' a is the number density of accelerated charged lep- 
tons and e e is the fraction of the kinetic energy of the falling 
baryons transferred to the accelerated leptons in the active 
region frame. 

In GRB/relativistic shock literature it is usually as- 
sumed that n' a ~ f'n' , i.e. only falling leptons are acceler- 
ated. However, the validity of this assumption is not certain 
because once the electric and magnetic fields are produced 
by the instabilities, all the charged leptons are accelerated. 
If the initial number density of charged leptons in the rest 
frame of the shells is similar to the slow shell, then as the 
flux of falling leptons is enhanced by a factor of 7' , the den- 
sity of local leptons can be neglected. However, in prompt 



4 We remind that column density and total power are scalars and 
therefore their value is frame independent 



collision one expects that the relative Lorentz factor of the 
shells be 0(1). In this case the local density of leptons is not 
negligible and n' a ~ {n' te + n' t p), with n' te and n' tp respectively 
the total number density of electrons and positrons. For a 
neutral matter with negligible positrons content n' a ~ n' te . 

Motivated by the power-law distribution of accelerated 
charged particles in other astronomical shocks, e.g. super- 
novae and cosmic rays, it is usually assumed that the dis- 
tribution of accelerated electrons responsible for the GRB 
prompt and afterglow emission is a power-law: 



Kile) 
N e 

7m 



7m 



-(P+1) 



for 7e ^ 7n 



2™ /2 

P m e n a 



n a p = 

7m (p — l)e e j' 2 m p n' 
(p - l)e e ~f' 2 m v n' 
pm e n' a 



(13) 
(14) 
(15) 



Recent simulations of particle acceleration by F ermi process 
in the relativistic shocks show (ISpitkovsky|2008T) that n' t ,( , y e ) 
is best fitted by a 2D Maxwellian distribution plus a power- 
law with an exponential cutoff: 



n'eile) = Ci7 e exp(-7 e /7i) + 



C 2 7e min[l,exp(-(7e-7i)/7e«t)] 



(16) 



where C2 — for y e less than a minimum value "/ m in- The 
typical values of parameters obtained from the fit to simula- 
tions for an initial Lorentz factor of 70 = 15 are: ymin = 40, 
71 = 6, 7, = 300, 7 cut = 100, and 5 = 2.5 (|Spitkovskvl 
2008). Implementation of this distribution makes the model 
presented here significantly more complex. In the range of 
energies relevant to the prompt and early afterglow emission 
of GRBs, the first term in (|16[l is negligible and for 7 e > 7; 
the distribution has the form of a power-law with exponen- 
tial cutoff. Conservation conditions similar to (|14|l and (|15|l 
for this distribution lead to the following relations: 



7m 

P 

Jcut 



-(p+1) 



min[l,exp(-(7 e -7i)/7 cut )] (17) 



1 - 


(7m) 







+ 



p( exp( — )r(-p, — 

7m / 7cut 7cut 



(18) 



-y' 2 m p n' e e 



N e j 2 



7m 



(P-1) 

O-i) 



1 - 



Ji_ 



-(p-1) 



+ 



exp (-lL)r(-(p-i),-^) 

"fcut 1 cut 



(19) 



where F(a, x) is the incomplete Gamma function. As the 
number of parameters in this distribution is larger than the 
number of conservation conditions, in contrast to the power- 
law distribution, it is not possible to find an expression for 
-/V e and 7 m with respect to the total density and the fraction 
of electric and magnetic energies transferred to leptons. A 
more simplified version of this model is a power-law with an 
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exponential cutoff: 
ne(7e) = N, 



7m 

for 7 e > 7, 



(p+i) 



exp(- 



7e 



7cui 
7cui ^> 7m 



(20) 



Note that due to the exponential cutoff the restriction to p ^ 
2 does not apply and the index of the power-law term can 
be negative. Using conservation conditions pip and (|12[l we 
can find relations between parameters of this distribution: 



7cui 



\p\m e n' a 

n a( — ) |r(-p, — ) 



(21) 
(22) 



7raJV e 

" 7m ' ' ' " ' 7" 

As we have only two constraints, it is not possible to find 
expressions for 3 constants N e , y cu t, and j m , and one of 
them will stay as free parameter. In Sec l2.4l we show that 
this type of electron distribution is necessary to explain the 
hard spectrum of the short hard and some of the long bursts. 

In the introduction we mentioned that the induced 
transverse magnetic field is produced by Weibel instabil- 
ity in the active region. The magnetic energy density is 
parametrized by assuming that it is proportional to the en- 
ergy density of in-falling particles to the shock front/active 
region: 



8tt 



e B c 7 m p n 



(23) 



It is expected that both e e and t B evolve with time/radius. 
If the central engine is magnetized the external magnetic 
energy can be very important and an external field should 
be added to the right hand side of (123[) . Here for simplicity 
we neglect such cases. 

Considering the simplest case of a power-law distribu- 
tion for electrons and also assuming that only in-falling elec- 
trons are accelerated i.e. n' a — n' y' (n' is the minimum of 
n' a for a radiative shock), the synchrotron term in © and 
((6| becomes: 



dE' 



4nmc 2 dr' 



4:aaTm,;n'oeie B y b r' 2 Ar' 
3m|/3' : 



p(p - 2) 

(24) 

For the reasons explained in the introduction and in Sec. 12. II 
we believe that in a realistic model of relativistic shocks one 
should consider the time evolution of electric and magnetic 
fields. Here we assume a simple power-law evolution with a 
constant index: 



(r'o) 



(25) 



Using (|24p . the expression for the column density (0, and 
the momentum conservation equation ©, we obtain the fol- 
lowing equation for the evolution of the relative Lorentz fac- 
tor: 



Nor'U 



1) + (3 - K )r' 2 n(r' )Ar'(r' W (l ~ %)P 



(3-«)(l-/?') 
4am 2 p a T N*e 2 e (r' )e B (r' W 6 r' 2 Ar' ( f 



The parameter rj = 2a e + a B + 2k is the evolution index of 



(26) 



the density and fields. Although this differential equation is 
of order one, it is highly non-linear. To solve it we proceed 
a perturbative method based on iteration. Moreover, it de- 
pends explicitly on Ar' and as we discussed in Sec l2.2l to be 
able to find an explicit solution for f'(r'), we have to model 
its evolution. We consider two models: 



2.2.1 Dynamically driven active region 

Assuming that the shock strength and consequently Ar de- 
pends mainly on the density difference, and that the densi- 
ties of the shells in their rest frame are roughly the same, we 
expect smaller Ar' for larger 7'. On the other hand when 
the relative Lorentz factor is small and the shock is soft, Ar' 
should be proportional to /3' and Ar' — > when /3' — > 0. The 
simplest parametrization of Ar' with this properties is: 



Ar' = An, 



7^ 



9(r' 



r ) 



(27) 



where Ar' is a thickness scale. A S-function is added to 
(|27|) to explicitly indicate that the expression is valid only 
for r' ^ r' . Note that in this model the initial thickness is 
not null and therefore it is assumed that it was formed in a 
negligible time or the value of Aro is the final value from a 
previous regime that makes the initial active region before 
(|27|) can be applied. This model should be suitable for the 
prompt/internal shocks in which two high density narrow 
shells pass through each other and one expect that roughly 
instantly a narrow and dens active region forms around the 
shock discontinuity (See also next section for other cases). 
As ft /"/ < 1 is expected to be a decreasing function of r', for 
t > the width of the active region decays and for r < 
it grows. But it is not always the case, see simulations in 
Paper II. 



2.2.2 Quasi-steady active region 

At the beginning of a strong shock, presumably an internal 
shock or when the slow shell is extended, roughly homoge- 
neous, and has a low density, we expect that after a transient 
time in which the active region grows, its thickness arrives 
to an stable state determined by the relative Lorentz factor, 
density, synchrotron emission, and expansion of the shells. 
This stability should persist until the loss of kinetic energy 
due to radiation and mass accumulation becomes important, 
or the fast shell passes through the slower one (this does not 
happen for a radiative shock). In this case we parametrize 
the time evolution of Ar' as: 



Ar' = Ar„ 



B(r - r' ) 



(28) 



where Aroo is the final width when the equilibrium is 
achieved. For the decay of the active region at the end of 
this regime one can use the dynamical model. Another pos- 
sibility is to consider: 



Ar D 



6(r 



r ) 



(29) 



In Appendix [A] we argue that with small modifications the 
calculation of dynamical evolution can be used for this model 
too. 
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This model is specially suitable for studying the exter- 
nal shocks of the ejecta from the central source with diffuse 
material or wind surrounding it, and/or the ISM. One ex- 
pects that in these cases the density of the relativistic ejecta 
- the fast shell - be much higher than the wind or ISM and 
its extension much smaller. In Paper II we show that such 
a model along with a late emission from internal shocks can 
explain the shallow regime observed in the X-ray light curve 
of the majority of GRBs detected by Swift . 

For studying the evolution of the active region, apriori 
we should also take into account the total size of the shells 
and the passage or coalescence time. However, observations 
show that synchrotron emission continues for a significant 
time after the end of the shock - when shells passed through 
each other or coalesced. For a far observer what is matter is 
the emission rather than physical encounter between shells. 
All these stages can be modelled by one of the models ex- 
plained here or similar models for the evolution of Ar'. In 
this case the difference between various stages of the colli- 
sion is reflected in the different value of parameters. 



2.3 Evolution of relative Lorentz factor 

To solve equation (|26[1 we use a perturbative/iterative 
method based on the assumption that the dimensionless 
coupling in the r.h.s of this equation is smaller than one. 



By dividing both sides of 
the coupling A: 



26)1 with nor one can extract 



d 



\ r ' i 



1 + O-.K^Ar'M) 7 > o(l _ W 



A 



(3-k)(1-0O 

2-1) 

l3>Ar>(r' ) \r'J 

4amla T Nl ) Ar'(r' )el(r' ) e B{r' ) 
3m'i 



Ai 6 Ar' (r 



(30) 
(31) 



It is straightforward to see that if the initial column density 
of the slow shell/ISM n' Ar'(r' ) < 10 22 cm" 2 , for any value 
of e 2 (»"o) < 1 and es(rg) < 1, the coupling A < 1. This up- 
per limit on the shell column density is in the upper range of 
the observed total Nh column density of GRBs. However, 
the real Nh can be much higher than what is measured 
from the absorption of the soft X-ray at least > 100 sec af- 
ter the trigger, because it i s in conflict with N h estimated 
from Lyman- a absorption (| Watson et aZ.|[2007T ) . The differ- 
ence can be due to the ionization of the neutral hydrogen by 
UV emission from the prompt emission. Nonetheless simu- 
lations of the formation of the electric and magnetic fields 
in the shocks show that the fraction of the kinetic energy 
transferred to the fields, specially to the magnetic field, is 
much less than one (|Keshet et aLl200Sr ) . Therefore even with 
larger column densities, the value of A should be less than 
one and the validity of the perturbative method is justified. 

The zero-order approximation corresponds to A — > 0. 
In this case eq. (|30[) is a pure differential and its solution is 
trivial: 



(3-k)T> 



-L O 



V = 



n'(r' )Ar\r' )(3' o7o 



N'r' 



(33) 



where (3'^ indicates the zero-order approximation for P'(r'). 
In the rest of this work we only concentrate on k / 3, but 
calculations can be easily extended to this exceptional case. 

The physical interpretation of (|32[) is quite evident. 
/3'(r') changes inversely proportional to the total mass of 
the shell including the accumulated mass of the swept ma- 
terial. This zero-order solution does not take into account 
the energy necessary to accelerate particles of the slow shell. 
Thus, its use without radiation corrections will lead to a vi- 
olation of energy conservation. Parameter T> presents the 
strength of the shock; smaller 23, faster the constant term in 
the denominator becomes negligible with respect to radial 
growth and /3( )(r') approaches a cubic decline (for « = 0) 
due to the adiabatic expansion. The origin of term (3 — k) 
is partially geometrical and partially related to the density 
variation with r in the slow shell. It is the effective mass 
accumulation index of the shock. 

As Q32p is the dominant component of the dynamics 
and is used through out this work, it is useful to have its 
asymptotic behaviour for (-p-) > 1 and (p-) 2> 1: 



E > 



(3-K)2?(^f) 
#(1- 



For (^) 3 - K » 



(3-«)P 



(3-k)U 



)(-H 3 

J\t> > 



(3-k)P 



> (^) 3 ~ K » i 

(34) 

We use the zero-order solution in the r.h.s. of (|30|l to ob- 
tain the first-order correction of the solution. The n-order 
approximation corresponds to using the (n — l)-order ap- 
proximation in the r.h.s. of (|30[) and solving the equation: 



v r'o ' 



1 + 



(3-K)n'(r-')Ar-'(r') 



^ 7o(l-/W(n 



(3-*)(l-0( n) ) 
A fiil^Ar^ 



Ar'{r' ) 



,2-77 



dx 



(35) 



(n-l) 



where T>r n ) is an integration constant. For any order of cor- 
rection (|35[) must satisfy the initial condition i.e. /3'( n \(r' ) = 
/3o and from this constraint one can determine 2?( n ). It is 
easy to see that ©( n ) = 2? f° r all orders of perturbation. 
Calling the integral term M.( n -i)(r') for (n — l)-order solu- 
tion, we find the following recursive expression for f3^ n y. 



Pi 



V 



AM, 



-iiO') 



n) 



_J_(f J±)3-K _\\ + 2r- AM 
3-k vV r '„ I ± ) T „, 



)(r') 



(36) 



At this point we have to consider a model for Ar'. For the 
dynamical model A1(o) becomes: 



M (0) (r') 



0(o 



Ar' 

..T-3+ 



dx 

0{ o) (.r') 



dy 
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The second form of the integral is obtained using ft^ solu- 
tion. Unfortunately this integral can not be determined ana- 
lytically. Nonetheless, by expanding one of the two terms in 
the integrand it is possible to find an approximation which 
is useful for getting an insight into the behaviour of ftm, its 
dependence on various parameters, and the importance of 
radiation correction. They are summarized in Appendix 1X1 
For quasi steady model Mio)(r') is: 



find following expressions for .M(o)(t~') when e <C 1 



M m (r') » ArV?(l--)(l + (3 --)/3 2 ) £ 
for \r/\ > 



(39) 



M m (r') = Ar' x I r '° X2 ~ V{ \. X ~% x - 



(l-^o) 2 ) 3 



3 — k 



((3- 



Jg' 



(0) 



dy 



y* 



rCi + (7 



(l-y*r y y (3- K )V ft 
l-((3-«)X»)-^(l+( ' 



(3 - k)V ft 



(38) 



The similarity of the second term in (|38|) to the integrand of 
(|36[) shows that at first order of radiative correction the two 
terms in (|28p act independently. The constant term leads 
to the first term in l|38p which is equivalent to (I36j> with 
t — > 0, i.e. an active region with constant thickness. The 
term proportional to (r' /r') s in (|28p is responsible for the 
second term in (|38p . Up to a constant it is also the expression 
for a constant thickness model with 77 — > 77 + 5. This means 
that the effect of power-law term in (|28[) is very similar to 
power- law dependence of fields and shell density on r'. This 
behaviour justifies the name quasi-steady we have given to 
this model. Analytical approximations of these integrals can 
be found in Appendix lAl 

In the last paragraph we considered power-law depen- 
dence on r' for the electric and magnetic fields. The case of 
an exponential rise or fall of the fields is important at the 
beginning and at the end of gamma-ray spikes or early X- 
ray flares. For this case, the only modification in (|37[l and 
(|38p is the replacement of x~ v with exp(— r/x) where here 77 
is a dimensionless coefficient determining the speed of expo- 
nential variation. It is negative for rising fields and positive 
when fields are declining, similar to the power-law case. In 
the same way if th e electron distribution n' K ('ye) includes an 
exponential cutoff (|Dempsev fc Duffvll2007l ') , an exponential 
term similar to the term for fields appears in the expression 
for Mm)(r'). Therefore in a general case, the integrand in 
()37[) and (|38p includes an exponential term which makes it 
even more complex. Nonetheless, the expansion of the expo- 
nential permits to obtain the analytical approximation given 
in the Appendix lAl 

Finally by using M(p){r') in (f36|) and ((9|, we can deter- 
mine /3(-n the first-order radiation corrected evolution of ft 
and column density n'(r')Ar' with radius/time. The com- 
plexity of expressions for M-(o)(r') and consequently for ft^ 
does not permit to investigate the effect of various quanti- 
ties from the exact calculations and we leave this for Paper 
II where we numerically evaluate the behaviour of shocks 
kinematic and radiation. Here we just consider the simplest 
cases when in (J34J), r'/r' > 1 or r'/r' > (3 - n)T>/ft . 
Using (|A6[) and ()A7[) respectively for large and small 77, we 



Ar'oB 



(3-k)(t-2 + 



3 — K 



+ 



A> 



, 2 /(3- K )(3-i)(r-l+3^) 



77(3-f)(l-l) 



7o 



(3 - n)ft 



for \r)\ > (40) 
(41) 

Po 



Therefore Mm(r') oc Ar'oBe for e < 1. The constant coef- 
ficient is expected to be of order 1. By applying these results 
to (1361) we find: 



#1 



V - AB'e 



W + 11 



AB')s 



(42) 



where B' is B multiplied by the corresponding constant coef- 
ficients in (I39p or Q40p depending on the value of r\. Compar- 
ing this result with the corresponding ft{ \ we conclude that 
the strength of the radiation correction of ft and its effect 

on the kinematic of the ejecta/jet depends on S = AB' . As 

— T? 

B is proportional to T> 3 -" , for a positive 77 and same A, 
larger T> (stronger shock), smaller <S. In this case the kinetic 
energy of the shock is much larger than radiation and there- 
fore the synchrotron emission does not significantly modify 
the evolution of the shock. Note also that although S is lin- 
early proportional to the synchrotron total coupling A, it 
depends non-linearly on T> through a power which depends 
on the time/radius variation of the electric and magnetic 
fields as well as the density of the shells. The quantity S in 
this model depends also on 70: larger 70, larger the influence 
of radiation. This simply means that the effective thickness 
of the active region decreases faster when the effect of ra- 
diation is stronger. This behaviour is consistent with the 
phenomenology of the model described here. 

Similar expressions can be found for the other extreme 
case i.e. when r' /Cr' ^> 1: 



. Ar' 07 oV 1 ^(l-^) 



M (0) (r') 



(3-*) a 



+ 



(3- f)A 2 , fCr' 3 -^ 



r + 1 + T 



+ 



(3-|)A 2 fCr' 3 ' K ^ 2 



for 1 77 1 > (43) 
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r + l 



3-n 



Cr 



(3-j)/T 



[1 + 



3—k 



+ 



(1 + 



(44) 
> 



^+1+3^ 

for |n > 

When r'/CrJ -> 0, M (0) (r') -> const., if r - 2 
for |n| > or r - 1 + ^ > for > 1. This can be 
interpreted as a saturation state for the synchrotron radi- 
ation which has been observed specially in bright bursts 
where peaks of the prompt gamma-ray emi ssion are roughly 
square -like. Examples are GRB 060105 dZiaeepour et al\ 
l2006ah . GRB 061007 dSchadv et al\ l2006bl . 120071). GRB 
0608 13 A jMoretti et aZ.ll200fJ ). and GRB 0704 27 dSato et al. 



2007alV Even the Super Burst GRB 080319B l|Racusin et al 
2008al lbh seems to consist of 3 overlapping square shape 



peaks. See also simulations in Paper II. 

The conditions mentioned above can be considered 
as consistency conditions because if they are not satisfied 
-M(0)( r ') — > oo which is not physically acceptable. There- 
fore these conditions constrain parameters of the model - 
rj, r, and k. For instance, for a slow shell with a roughly 
constant density k ~ 0. Therefore n depends only on the 
behaviour of electric and magnetic fields. If the variation in- 
dex of the these fields are small and positive, the value of r 
can be small, i.e. the radiation will not vary very quickly. By 
contrast, a negative index - increasing fields - can not last 
for long time and imposes a large value for r. This simple 
argument shows that there is an intrinsic relation between 
these parameters. However, only a detailed modelling of the 
microphysics of the shock will be able to determine possible 
relations and their physical origin. Nonetheless the discus- 
sion above shows that the simple model studied here is con- 
sistent and includes some of the important properties of the 
phenomena involved in the production of GRBs. 

Using (|36[) one can see that a constant M (o) (r') is equiv- 
alent to redefinition of T>. Therefore when the radiation term 
arrives to its maximum, evolution becomes like a non- 
radiating ejecta. 

In the case of a quasi-steady active region, the behaviour 
of M(o)(r') is essentially similar because of similarity be- 
tween two models explained above. However, the time scales 
and indexes are different. For instance, when e <C 1 the value 
of -M(o)(r') is proportional to e with a coefficient equal to 
(|39[) or (|40[) and r = 0, minus the same term with r; — > rj + 5. 
In this case, the initial y' does not have an explicit contri- 
bution and the slope of _M(o)(V') is smaller than dynamical 
model. For r' /Cr' 3> 1 if other parameters are the same as 
dynamical model, the absence of 7q and smaller power of 
Cr' /r' means that .M(o)(r') approaches its maximum value 
slower. There is also a slower change when one of the two 
contributor terms becomes too small and negligible. As ex- 
pected, all these properties make this model more suitable 
for modelling the afterglow. 



Up to now we have only discussed the solutions of dy- 
namical equation (|26[1 corresponding to the rise of the syn- 
chrotron emission. The falling edge of the emission i.e. when 
^(0)( r ') ~* can be obtained simply be time/radius re- 
versal of rising solutions, see eq. (|29p . For instance in (|39|) 
and (|40p if r' < r' , e < 0. By moving the initial condition 
from r' to r', we obtain a positive but decreasing value for 
A4(o)(r') which becomes zero at r' — r' . These approxima- 
tions however do not permit to determine when the radiation 
begins to decrease. For this we need a detailed study of the 
evolution of fields and other shock properties. 

In summary, the evolution of f3' determines the kine- 
matic of the burst and is important for the estimation of all 
observables such as what we will discuss in the next sections 
- synchrotron flux, hardness ratios, etc. is also impor- 
tant for determining the evolution of other parameters that 
are not directly observable and a model must be used for 
their extraction from data. A good example is the time vari- 
ation of uj' m the minimum characteristic frequency of the 
synchrotron emission. It determines the behaviour of the 
spectrum and light curves (see expression (|78[) below for its 
definition). Assuming the simplest case of n' a — n' in (|15p . 
from the definition of u>' m one can see that uj' m oc j' 2 eB- The 
proportionality coefficient is time/radius independent. 



2.4 Synchrotron flux and spectrum 

In this section we first remind the synchrotron emission for 
the purpose of completeness and then we use the results for 
determination of lags between the light curves of different 
energy bands. 

The ejecta from a central engine that produces the 
gamma-ray burst is most probably collimated and jet-like, 
otherwise the observed energy is not explainable. On the 
other hand for the far observers, even a spherical relativis- 
tic emission looks collimated to an angle < 1/F along the 
line of sight where V is the Loren tz factor o f the emitting 
matter in the observer rest frame (|Reealll967h . Therefore in 
any case we need to consider the angular dependence of the 
synchrotron emission. Moreover, we need to consider the de- 
lay as well as angular dependence of the Doppler shift of the 
emission. Simulations show that even with an angular inde- 
pendent emission, these effects can apriori explain the lag 
betwee n different bands observed in both BATSE and Swift 
bursts (|Qinll2002l ; TOm et aZjIaxii ; iLu et fflHl2006T l. However, 
in these simulations the spectrum and the time profile of 
emission have been put by hand. 

There are a number of evidence against a high lati- 
tude/Doppler effect origin of the observed lags. First of all 
the total effect of high latitude emission decreases when 
Lorentz factor is very high. It is expected that in GRBs 
r > 100, and therefore this effect should be very small. In 
addition, even in the early X-ray emission where it was ex- 
pected that this effect dominates, it has not been observed. 
On the other hand, for a given category of bursts, short or 
long, it does not seem that there is any relation between lags 
and spectrum as expected from a Doppler effect. Therefore 
we conclude that the contribution of high latitude emission 
and Doppler effect is sub-dominant and can not explain the 
observations. On the other hand, we show that even with 
neglecting the high latitude emission, due to the evolution 
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of physical properties a lag between different energy bands 
exists. 

In the model presented here the synchrotron emitting 
matter is confined to the active region. Therefore we identify 
the bulk Lorentz factor T with respect to the observer with 
the average Lorentz factor of the active region with respect 
to the observer. It can be related to the relative Lorentz 
factor 7 '(r') obtained in Sec. l2.2l and the final Lorentz factor 
when two shells are coalesced: 



r(r)=r /7 '(r)(l + /3 / /3'(r)) 



(45) 



where F/ is the Lorentz factor of coalesced shells with re- 
spect to observer. The initial value of 7 '(r) is the relative 
Lorentz factor of the colliding shells. In the case of an ex- 
ternal shocks on a low velocity surrounding material or the 
ISM Tf ~ 1 and T(r) « Note that here we have writ- 

ten 7 with respect to observer coordinate r because on what 
concerns the synchrotron emission, only the observations of 
far observers matter. 

The e nergy (intensit y) angular spectrum of synchrotron 
emission (|jacksonl l200lh by one electron or positron in a 
frame where it is accelerated to a Lorentz factor of 7e iijfj 



d 2 r 

dw'dQ' 



41T 2 LU' C C 

'2a' 2 



7e 



* 1/3(0 



c 



1 + 7 I6 
3 3 /c\ = 3ej 2 e B' 
2 Je \p'J 2cm e 



(46) 
(47) 
(48) 



Quantities with a prime in (|46|l to (|48f) are with respect to 
the frame where Lorentz factor of electrons is 7 e . Here we 
identify this frame as the rest frame of the active region, p' 
is the Larmor radius of electrons. The angle 9' is the angle 
between acceleration direction and emission. Without loss 
of generality it c an be assumed to be in x-z plane. Therefore 
dO. = cos9d9d(f> (|jacksonll200ll '). When \9\ > 0, dO. w d6d(j>. 

To obtain the power spectrum that is the measured 
quantity, we must divide the intensity by the precession pe- 
riod of electrons 2irp' /c. We expect that accelerated elec- 
trons have a range of Lorentz factors, therefor we should 
also integrate over their distribution: 



d 2 P' 



Lo'dio'dQ' 



4 ^3 C / d > <(7e)7e 1 f ^7 ) (1 + lid' 2 ? 



K | /3 (C) + 



..2 a/ 2 



i + ii 



S72 K 1/3(0 



(49) 



Note that we have divided the angular power spectrum 
by u' to make it dimensionless. I n the observer's rest 
frame the spectrum is transferred as (|Sadun fc Sadunlll99ll : 



5 The expression l|46|l is valid for small 9' angles. However, the 
main part of the emission is in "y e 6' < 1, and the intensity expo- 
nentially decreases fo r ■y^d 1 > 1. Th erefore it is a good approxi- 
mation for any ang lc (Jackson 2001). 



iHuang et qi1l2000l ; Ivan Eerten fc Wiiers|[2009l ): 

d 2 P 1 d 2 P' 



LudtodQ 



r 2 (l - /3cos6>) 2 uj'duj'dQ,' 
t 



cos 9 - 



r(l -/3cos( 

cos 9' + p 
1 + /3cos 9' 



cos 9' = 



cos 9 - (3 
1 - j3 cos 9 



(50) 
(51) 

(52) 



where in equations (|50|l to (|52|l . (3 is related to F the Lorentz 
factor of the active region with respect to the observer. To 
find the total power at a given frequency P(t, uj) in the frame 
of a far observer, we integrate over the distribution of accel- 
erated electrons and the emitting volume but constrain it to 
the emission in the direction of the observer. As we assumed 
that the active region is thin, we neglect the absorption of 
synchrotron photons inside the active region itself. Without 
loss of generality we put the observer at f2 = Q! = 0. In 
this case for an electron moving at angle fi^ with respect to 
observer, only photons emitted in the direction of Q' — Q[ 
are detected by the observer. Therefore we need to integrate 
either on fi^ (or equivalently Qi) or on il'(Q,). For simplic- 
ity of notation we use the latter. As the synchrotron angular 
distribution does not depend on 4> ( or 0') we only need to 
integrate over 9 (or equivalently 9'): 



dP 

ujduj 



2ir 



dr 



d9 cos 9 



d 2 P(t- At),9,w) 
uidudQ, 



(53) 

9min and 9 max are minimum and maximum visible angle for 
the observer with the constraint \9 m i n \ < 1/T and \9 max \ < 
1/r. We define A9 = \9 rna x + 9 m i n \/2 as the view angle of the 
observer with respect to the ejecta/jet axis. These angles are 
not directly measurable and therefore the simplest assump- 
tion is a symmetric ejecta 9 max — —9 m i n = 1/r, i.e. A9 = 0. 
Note that photons coming from 9 ^ arrive to the observer 
with a time delay At where At(9) is (|Rvbicki fc Lightmanl 
l2004h : 



At(9) = 



r(l — cos ( 



(54) 



c/3(r)(l + /3 2 F 2 sin 2 9)2 

In (|54p we have assumed that the initial radius from which 
the shells are ejected from the central engine is much smaller 
than their distance from it when they collide. Therefore, the 
initial radius is neglected. Due to the direct relation between 
radius and time in this model, t — At(9) can be replaced by 
Ar(9) = r- c(3(r)At(9). The quantity Ar(9) should not be 
confused with the thickness of the active region Ar. 

In this model we have considered Ar r and the physi- 
cal properties of the active region are close to uniform. Thus, 
the integral over the interval r and r + Ar becomes trivial if 
we consider r to be the average distance of the active region. 
This is similar to the way we calculated kinematic quanti- 
ties in Sec l2.2l With this simplification the total spectrum 
becomes: 



dP 

iudui 



2nAr 



d9 cos 



a d 2 P(r- Ar{9),9, 
cudujdQ 



(55) 



If we use (0 to describe r as a function of time, equation 
()55|) depends only on t and u>. 

The width Ar can in general depend on the energy. In 
Sec l2.ll we described that when electrons lose their energy, 
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they get distance from the shock front or in another word 
are pushed to the upstream. Although in this region the 
magnetic field is expected to be weaker, it can be enough 
for the low energy emissions, UV, optical, IR. Moreover, in 
a collimated or structured je10 these less accelerated par- 
ticles emit mostly in lower energies. In this simple model 
of shock one way of taking into these effects is to consider 
that Ar as well as 8min and 9max depend on the energy. For 
Ar we can simply assume that Aro in dynamical model or 
Aroo in the quasi-steady model are energy dependent. Es- 
timation/modeling of energy dependence of 8 m i n and 8 m ax 
is more difficult because the only way to modify them is 
through r. Assuming 8j e t > 1/F we must consider that V is 
^-dependent. This needs a modification of the dynamics and 
makes the model too complicated. For this reason here we 
ignore the energy dependence of opening angle of the jet. 

The differential term in (|55[l can be replaced by (|50p 
and (|49[) . We must also take into account the proper time 
delay discussed above. Therefore: 



the total energy spectrum can be written as: 



dP 

uiduj 



r 2 Ar( — 

2n* T Wcc/Je r . 

r(l — cos 9) 
(l+/? 2 rsin 2 0)1/2' 



d6 cos 9 



r 2 (r 



(l-/?(t) cosOf 



■(l+7eV 2 ) 2 



2 a a 



Kl 3 (C) + 



7 c 



1+7, 
lo' 3eB' 



ft 



2cm e 



(56) 
(57) 



As the angle 9 and therefore the delay is small, we use Taylor 
expansion around r to obtain the explicit expression for F 
and j3 in the integrand of (1561) : 



r( r(i-cosfl) ) « r(r)(i - 

l (l + /3 2 r 2 sin 2 0)V 2J 

r(l-cosfl) rfr(r) 



+ ...) 



(1 + /3 2 r 2 sin 2 (9) 1 / 2 Vdr 
In a similar way we can expand (3(r, 8) around 
r(l — cos C ' 



(58) 



^ ( (IT^w^ )w/?(r)(1 " 

r(l-cos0) d/3(r) 
(1 + /3 2 F 2 sin 2 6*) 1 / 2 I3dr 



+ ■ 



(59) 



As | sin(?| ~ 9 ^ 1/r, we can also expand the 9 dependent 
terms in (|58[) and (|59|l . It is more convenient to use the 
relation between 9 and 9' in (|52p and transfer 8 to 8' . We 
keep only terms up to 8' 2 order. With these simplifications 



B A structured j et is usually considered to have a tr ansverse gra- 
dient in density iXu et aZ.1 '2005; Tak ami et "alll2007t) . However it 
is expected that gradually a transverse gradient in Lorcntz factor 
forms too. 
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P' l3'(r) 
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(60) 



(61) 



The coefficient Q(r) in (|60|l presents the lowest order cor- 
rection due to the Doppler effect and delay of high latitude 
emissions. As expected, for V — > oo, G(r) cx 1/r 2 — > 0. More- 
over, G(r) does not depend on ui' (or equivalently to), and 
therefore it does not affect the lag between different energy 
bands itself. 

Apart from the energy dependence of Ar, 9' min , and 
9' max mentioned before and has an intrinsic origin, i.e. are 
related to the evolution of the active region and the structure 
of the jet/ejecta, there are two other effects that can create a 
lag: angular dependence of the emission and time/radius de- 
pendence of 7 m and ui' cc . Assuming a homogeneous density 
for the colliding shells, equation (|15|) shows that the time 
dependence of 7 m is mainly due to variation of the relative 
Lorentz factor 7' and the electric field energy fraction in- 
dex e e with time. The variation of ui' cc is also possible if the 
magnetic field changes with time. If fields are considered to 
be constant during the emission - as it is the case for many 
GRB models in the literature - only dynamical friction and 
decrease in 7 remain. They may be insufficient to explain 
the observed lags because the Doppler and high latitude 
correction is usually small, specially when T is large. In fact 
simulations show that with an angular independent emission 
flux, large lags can be obtained only if the emission happens 
at l arge radius wit h respect to the central engine, > 10 15 
cm (|Lu et al.ll20"0lih . This is orders of magnitude larger than 
what is expected from internal shocks. As synchrotron emis- 
sion is highly angular dependent, lags are expected not only 
from Doppler and high latitude correction but also when 
these effects are ignored. In fact this can be seen in (|60p . 
The integration over 7 e can be expressed as /(Cm) where 
Cm is C (defined in l|48|)) for "/ e — "/ m . It is then clear that 
r and uo' dependence of the integrand are not factorisable, 
and therefore even when Q(r) is ignored, the r dependence 
of dP/uduo can not be factorized. Therefore, giving the ex- 
pectation of high collimation and large Lorentz factor, the 
main cause of the observed lags in the GRBs seems to be 
the time variation of physical properties of the active region 
such as Ar, j m , and electric and magnetic fields. 

We note that the zero-order term in H60j> is the same as 
the expression (|49[) for the synchrotron emission by one par- 
ticle. Therefore, for this part of integration i n (1601) we can use 
the well known e xpression for t h e spectrum j Schwingcr 1949; 
IWestfield|[l959l : iJacksonl l200ll : iRvbicki fc Lightmanl 120041 ). 
For determining the Doppler and high latitude correction 
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term we can again use the integration method used for the 
zero-o rder term and express (|60|l as a sum of Bessel func- 
tions [|Westfieldll959h . Finally we find the power spectrum 
with the first order correction for Doppler effect and high 
latitude delays^ 



dP 

uiduj 



V3e 2 2 Ar 

T 

3tt r( 



^ / d7 e n' e (7 e )7e 



22 



T #i/s(-r)-9 — # a /s(— ) + 



(62) 



We can further integrate (|62[) if we consider a specific distri- 
bution for accelerated electrons n' e ( , y e ). In Appendix [Bl we 
calculate dP/uidu} for the power-law distribution (|I3p and 
comment the case for a power-law with an exponential cutoff 
distribution. Here we use those results to discuss the lowest 
order properties of the spectrum and light curves. 



2.5 Lag between light curves 

The presence of a lag between variou s energy bands of 
the G RBs has been det ected in BATSE dKouveliotou et all 
QUI) and INTEGRAL (|Folev et all 12003) light curves and 
confirmed by Swift -BAT (|Sakamoto et all 12007). Observa- 
tions show that in most long busts there is a significant lag - 
from tens to few hundreds of milliseconds between soft and 
hard bands. By contrast, in short bursts the lags is very 
small and within the present sensitivity and time resolu- 
tion of gamma-ray telescopes it is consistent with zero or at 
most a few milliseconds. These two classes are respectively 
associated with the explosion of massive stars (collapsars, 
hypernova) and to the merging of two compact objects - a 
neutron star with a black hole, two neutron stars, or a neu- 
tron star and a white dwarf. There is however evidence for 
the difference between lags of the separate peaks in the same 
burst. Moreover, some bursts that accordin g to their T90 can 
be cla ssified as short such as GRB 080426 ( Ziaeepour et all 
2008b) have relatively long lags of few tens of milliseconds . 
On the other hand some appa rently long bursts s uch as GRB 
0606 14 l|Mangano et alj|20u7l ) and GRB 080503 (|Mao et all 
2008) have small lags similar to the short bursts. Various 
explanation have been put forward for these out of norm 
behaviours: sensitivity of detractors only to the peak of a 
long burst leading to its misclassification as short, existence 
of a separate class of GRBs with intermediate durations and 
lags, long tail emission in otherwise short burst for bursts 
with long Tgo and small lags. Some authors even rule out the 
association of long GRBs - hypernova, short GRBs - colli- 
sion of compact objects, and suggest that they should be 
classified according to their lags: Short lags old population, 
long lags young population l|Zhang et al . 2007a|). Apart from 



7 For simplifying the integration over 8' in (1606 . we consider 
cos 9' + /? Si 2. This is a valid approximation when /3 — ► 1 and 



classification of progenitors of the GRBs, lags along with lu- 
minosity have been also used as pr oxy for the GRBs redshift 
determination dNorris et al.ll2000l ). 

In summary, lags seem to be important quantities re- 
lated to the nature of central engine of the GRBs, proper- 
ties of the ejecta, and the surrounding material. They can 
be relatively easily measured, and therefore it is important 
to be able to relate them to these phenomena. Some au- 
thors have tried to explain lags just as a geometrical effect 
relate d to the h i gh latitude emission and associ ated Doppler 
shift (|Qinll2002l ; [Qin et aZ.|[20o3 : iLu et ffllll2006h . As we have 
discussed in the previous section and also regarding r and u 
dependence of the spectrum in (|62jl . it is evident that even 
without Doppler shift and high latitude corrections the r 
dependence - equivalent to time dependence in the model 
discussed here - and u> dependence are not factorisable, and 
therefore light curves in different energy band can not be the 
same even when they are normalized to smear the amplitude 
difference. 

To define and determine lags we need fast varying fea- 
tures such as peaks. In the gamma-ray bursts light curves 
peaks are mostly observed in the prompt gamma-ray energy 
bands. Nonetheless, fast slew of the Swift satellite and some 
of the ground based robotic telescopes have permitted to 
observe the counterpart peaks or at least evidence of their 
presence in X-ray and optical bands. A realistic model for 
the lags should be able to predict the lag in all these energy 
bands if they have the same origin. On the other hand, a 
deviation of some bands from predictions can be used as an 
evidence for a different origin of the corresponding feature. 

The nature of peaks and their profile are not well under- 
stood. In the framework of the synchrotron emission from 
internal shock model as the origin of the prompt gamma-ray 
emission, the rising side of a peak indicates the beginning 
of the collision between shells and formation of the electric 
and magnetic fields that leads to the acceleration of charged 
particles and synchrotron emission in the induced magnetic 
field. Decreasing edge corresponds to the separation and/or 
the total coalescence of the shells. However, it is expected 
that even before separation/coalescence, microphysics in the 
active region arrive to a roughly steady state during which 
only slight changes due to e.g. density fluctuation in the 
shells will occur. In particular when the initial evolution of 
the microphysics is much faster than the time of the passage 
of the shells through each other we expect that for a limited 
duration the active region has quasi steady characteristics. 
Assuming such a case - in accordance with the discussion 
about the evolution of Ar in Sec. 12.21 - a peak corresponds 
to: 



d ( dP 

dt \ ujduj 



c/3(r) 



d ( dP 



dr V ujdu) 



(63) 



i dominates the synchrotron emission. 



The lag for a given peak corresponds to the difference be- 
tween peak time/radius for two frequencies or two energy 
bands. Usually observations are performed in known energy 
bands. Therefore, the purpose of the lag measurement is to 
determine the difference between r pea k the solution of equa- 
tion (|63[) . at two different energies. In the rest of this section 
we use the results of Sec. !2.4l to determine the lags. 

Using (|62l) the peak equation (|63[) can be written as a 
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partial differential equation: 
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(64) 
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(66) 



Note that the functions D and Q' depend only on r and not 
on uj' . If Ar is energy independent, so is IC(r). The spectrum 
in (1621) can be written as: 
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With this definition the partial differentials of dP/ujdui in 
the peak condition equation (|64[) can be calculated: 
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(71) 
(72) 



This equation is obviously very complex and solving it is 
not a trivial task. However, we are only interested in the 
change in the roots with respect to energy uj' . Moreover, if 
we restrict ourselves to \Auj' /uj' \ = \oj[/uj' — 1| < 1 and 
assume that in this case the corresponding roots ro and r\ 
are close i.e. [ri/ro — 1| < 1, then we can expand functions in 
(I72[) around uj' and determine the lag in the observer frame 



i.e. cf3At = n ~ ro. Assuming that w 1 we find: 

oujor \ Loduj ' am' or v ujcLuj / 



cAt 



a 2 i dp \ 

dr z V ijjduj ' 



P(r ,u' ) 



d 2 / dP \ 
dr 2 \ i^jdui ' 

Auj' uj'P(r ,m' ) 
lo' Q{r ,uj' ) 

K(r ) J d7 e nU7e)7e" 2 ^ - 

d 

D ( r o)-^-;(n'e('Y™h™ 2 'H) +Q'(ro) 

_ 2 (>')l 1 d I <\P 



(73) 



d7e^'e(7e)7e » 75^ 3~ ~7~ 

ouj'oQ uj' T dr \ ujduj 



(74) 



Q(ro,Wo) 



-^( r o) J «4(7e)7e 2 H + £(r ) 



^'(ro) 



W e(7s)7e 2 ^ - £>ne(7m)7m 2 ?i 



as 



-^-{Dn' e ('y m )^ m 2 H) - 5'(r )I>n / e (7 m )7 m 2 ^ + 



dg' , 

dr 



(r ) / n' e (7e)7; 



i gTj JC d / dP 
dg T dr\ uiduj 



(75) 



All the terms in (|74|l and (|75[l can be expressed as a sum 
of Bessel functions and their integrals using the definition 
of TL. For the special case of a power-law distribution of 
electrons most of the integrals can be calculated analyti- 
cally. We present the results in Appendix [B] Note also that 
Aw'/w' = Auj/uj and uj'd/du' = ud/duj. Other uj dependent 
terms also are F-independent and therefore the calculation 
of the lag does not need a pre-knowledge of the bulk Lorentz 
factor r. 

The left hand side of (|63[1 or equivalently (I64p is the 
slope of the spectrum (light curve) and therefore its en- 
ergy dependence explains for instance why peaks or more 
preci sely their auto-corre lations are wider at lower ener- 
gies ()Fenimore et aZ.lll995l ) . However, the complexity of ex- 
pressions (|73[) to (|75|) make the analytical estimation of the 
time and energy dependence of light curves and lags diffi- 
cult. In Paper II we present some simulations in which the 
lags are consistent wit h the long bursts or a re very small 
similar to short bursts (Saka moto et all 2006). 



2.6 Break at low energies 

Before finishing this section we want to make a few com- 
ments about the jet break which is one of the most impor- 
tant aspects of the afterglow light curves predicted since the 
beginning of the modelling of GRBs. 

Since the early days of the discovery of GRBs and mea- 
surement of apparently huge amount of energy released in 
these phenomena - 0(0.1 — 10) x 10 53 ergs for long bursts 
or even larger in a few except ionally bright bursts such as 
GRB 990123 GRB 080319B dRacusin eV aZ.I l2008al lbh and 
GRB 080607 (Ma ngano et al . 2008) - it has been suggested 
that these measurements are biased by the collimation of 
the fireball and the actual total emitted energy should be 
much smaller. To produce the gamma-ray prompt emission 
in a shock the ejecta should be highly relativistic with a 
bulk Lorentz factor of order of few hundreds. This leads to 
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a strong apparent collimation of the radiation from a spher- 
ically symme tric e.jecta to angles O < Oboost = 1/T for a 
far observer l)ReesHl967h . It is also possi ble that the e.ject a 
is not intrinsically spherical but jet like (|Rees et al.l ll998). 
In this case after deceleration of the fireball during its prop- 
agation through the surrounding material or the ISM, at 
some radius/time Qboost > open = 9 max — min and ra- 
diation is n o longer collimat ed. This lead to a drop of ob- 
served flux (|Sari et aZ.1 1 19991 ^1 . As this effect is purely kine- 
matic/geometric it should not depend on energy. Observa- 
tions of the Swift and robotic ground telescopes however 
contradict this expectation. The breaks seen in the afterglow 
light curves are usually chromatic. In many bursts no opti- 
cal break has been ob served up to million s of seconds after 
the prompt emission <|Roming et al\ fo008). In some bursts, 
mostly bright ones but not always, no br eak has been ob- 
served in X-ray either (|Sato et oZ.ll2007bh . In a significant 
fraction of bursts multiple breaks in the X-ray light curve 
have been observed. Apriori only one of these breaks (if any) 
can be due to the jet break. Therefore, the mechanism of 
break is more complex than just a kinematical effect. Here 
we want to argue that even in the simple case of the opening 
due to deceleration, we should not expect to have an achro- 
matic break if the afterglow is a synchrotron radiation. 

For obtaining the analytical expression (|62[1 for the ra- 
diation spectrum, we made the simplifying approximation 
that |0 m in| = |^maa;| — > oo . This permitted to analytically 
integrate the angular integral in (|60|l . Without this approxi- 
mation (|60[) depends on the opening angle. The justification 
for this approximation was that the synchrotron emission 
is highly directional and for 8'"/ e 3> 1 the intensity reduces 
exponentially because the value of £ in the integrand of (|60|) 
grows. £ increases with energy and therefore high energy 
photons are preferentially emitted at small angles with re- 
spect to electrons boost direction. When the collimation is 
reduced, the observer re ceives less high energy photons from 
high latitude electrons (|Panaitescu fc Meszarosl 1 19981 ) and 
therefore the effect of the jet break influences high energy 
light curves earlier than low energy ones. Conceptually this 
effect is very similar to the lag between light curves. It is 
however more difficult to investigate it mathematically be- 
cause the angular integral (|60p can not be determined ana- 
lytically as explained above. 

On the other hand, it is not sure that the reduc- 
tion of Lorentz factor can explain the absence of a break 
in optical frequencies when it is observed in X-ray. A 
number of su ggestions have b e en put forward to explain 
observations (|Panaitescul 120051 ; iMisra et all 120071 1 . In the 
frame work of the model presented here it can be at least 
partly explained by the energy dependence of the geom- 
etry of the active region. In fact the material that has 
lost its energy and has been decelerated is pushed behind 
the front edge of the shock. Due to the scattering it also 
has a rel atively larger transv erse momentum (See simula- 
tions in (|Vigelius et al\ I2007T) for non-relativistic shocks). 
Thus, less boosted particles in the ejecta have intrinsically a 
larger opening angle. This leads to a later jet break for the 
softer photons originating in this extended region. Therefore 
one does not need add i tional processes dRees fc Meszarosl 
19981; ISari fc Meszarosl |2000D. energy ilDai fc Lul 1 19981 ; 



The complex geometry of matter and fields in the mate- 
rial and shock can explain apparently contradictory obser- 
vations. As for multiple breaks, they can be due to the fact 
that multiple shells are ejected by the central en gine. This 



is in the same spirit as t he structured jet models ( Xu et ali 
l2005l ; lTakarm~ et aZ.| [2007). If they do not completely coalesce 
both the tail emission and the afterglow due to the external 
shocks will combinations of radiation from separate remnant 
shells. In this case the break of radiation from each shell is 
independent of others and a far observer who detects the 
total emission will observe multiple breaks in the emission. 



3 EXTRACTION OF PARAMETERS FROM 
DATA 

The model described in the previous sections has a large 
number of parameters. It would not help to better under- 
stand the physics of gamma-ray bursts production and their 
engine if we can not estimate the parameters of this or any 
other model from data. On the other hand, the only con- 
veyor of information for us is the emission in different energy 
bands. Therefore, we must find the best ways to extract the 
information as efficient as possible. In this section we sug- 
gest a procedure for extracting the parameters of the model. 
We also discuss their degeneracies. 

Since the massive detection of GRBs by BATSE, many 
efforts have been concentrated to understand their spec- 
trum. However, very little information could be extracted 
from the spectra. The complexity of time and energy depen- 
dence of the spectrum (|62[) explains why the time averaged 
spectrum does not carry extractable information. As GRBs 
evolve very quickly both in time and in energy, integration 
over these quantities smears the useful information. Unfortu- 
nately the effective detection surface of present gamma-ray 
detectors is too small and we do not have enough photons 
to make a broad spectrum in small time intervals. In this 
situation hardness ratios and their time variation are more 
useful quantities. Therefore we first discuss what we can ex- 
tract from hardness ratios. 

Consider a power-law distribution for electrons, equa- 
tion (|B1|) shows that at a given time the spectrum can be 
expanded as a power of u>' /tJc C 7m- The coefficient in front 
of the integral term does not depend on energy. Therefore, 
if we neglect Q(r) i.e. assuming the bulk Lorentz factor is 
very large, the hardness ratios depend only on the integral 
determined in (|B4|I . We define the hardness ratio between 
two bands with logarithmic mean energies U\ and u>2 as: 
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HR 12 (r) = — « ^ x ^ (76) 



Rees fc Meszarosl l200d: IZhangfc Meszarosl l2002h and /or 
components ( De Pasauale et a l. 2008) to explain the breaks. 



A max min tr7t7\ 

0Ji=<jJi - ui t (77) 

The approximations for integrals are valid when Auji/uji <C 
1. Expression ()B4|) shows that at a constant time/radius 
hardness ratio depends only on p, electron distribution in- 
dex, and on the Lorentz invariant factor u' /u)' m where: 

win = <4c7m (78) 

is the minimum synchrotron characteristic frequency. There- 
fore by fitting observational data with ([76} one can obtain 
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the index of the electron distribution and the time variation 
of cu m which is a very important quantity. Note that as HR 
is time/radius dependent we need multiple energy bands and 
their hardness ratios to make a statistically meaningful fit. 
If we release the assumption of power-law distribution for 
electrons, we can use (|62p to get insight into the distribu- 
tion of electrons Lorentz factor and ui m , but they will be 
degenerate. 

Once the energy and temporal behaviour of the integral 
term in (|62[) is found, the extraction of leading time depen- 
dent coefficient T(r) (defined in (J68J) ) from light curves is 
easy. Then if we use one of the models considered in this 
work or another model for the evolution of Ar, we can de- 
termine time/radius variation of Y{r) the total Lorentz fac- 
tor of the active region. Using (|45|l and (|36|) one can apriori 
estimate M(r'). However, the latter depends on the relative 
Lorentz factor j'(r'). If we can roughly estimate the end of 
the collision from the form of the light curve - for instance 
if we assume that the end of the collision is when the con- 
tinuous component begins an exponential decay - then: 



I> = T(r) 



(79) 



where 77 is the radius at which the coalescence of the shells 
finishes. 

w'/i, A4(r'), and 7 m (r') depend on the fundamental 
quantities such as magnetic field and thereby on the frac- 
tion of kinetic energy of the shell transferred to the electric 
and magnetic fields, the density of shells, and the initial 
distance of the collision from the central engine. But with 
only 3 quantities mentioned above we can not determine 
all these quantities and their time/radius variation. An im- 
portant observable that can help to extract more informa- 
tion from data is the lag. If we neglect the term propor- 
tional to G{r) and its derivatives, we can determine functions 
P(r ,uj') and Q(r 0l cu') in and (J75J from the knowl- 

edge about T(r) and the integral of H(r) (defined in ((U^J]l ) 
explained above. Thus in this way we can determine the 
lag between two energy bands. Comparison of the observed 
lags with the expectations from the model helps to estimate 
some of other quantities such as tq. However, as the ana- 
lytical expression for the lags is very complex, only simula- 
tions and fit can permit to solve this inverse problem. Using 
a simpler version of the model presented here, this proce- 
dure has been applied to analy se the Swift data for GRB 
060607a dZiaeepour et fflZ.ll2006bh and to explain some of its 
peculiar properties ( Ziaeepour et al.l l2008a). 

At the beginning of this section we mentioned that the 
most complete information about GRBs is in the spectrum 
at short time intervals when the flux does not change signif- 
icantly. In fact in practice hardness ratios are calculated by 
adding together photons in a given energy band and time- 
rebinned event data to reduce the noise. When multiple si- 
multaneous hardness ratios are available, they can be con- 
sidered as a low resolution normalized spectrum in a short 
time interval. Because these quantities give the most direct 
insight into the physics of the collision, in the rest of this 
section we investigate in more details the spectral behaviour 
of the flux at constant time. 

First we consider a power-law distribution for electrons. 
When Lu' m > u' hypergeometric functions 1F2 in ()B4[1 and 
(|B5[) can be expanded as a polynomial with positive or zero 



power of ui'/uj' m . This shows that at the lowest order in 
ui' /ui' m the spectrum dP/ojduj oc (a/ /io' m )~ 2 ^ and therefore 
HR12 oc (u^luiT 2ri . The power of the coefficient in the 
second 1F2 term is larger and therefore for determining the 
spectrum at zero order it can be neglected. But for the first 
order expansion of (|B1|) this term is significant and when it 
is added the spectrum becomes: 
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(80) 

The amount of deviation from the dominant power-law only 
marginally depends on the index of the electron distribution 
function p. When uj'/u)' m — > 1 the higher power of u>' /u}' m be- 
come important and should be considered. However, as they 
are all positive, they make the spectrum harder. This re- 
sult shows that although small hardening of the spectrum 
can be obtained in this regime, it is not possible to have a 
softer spectrum with an index smaller than —2/3 and closer 
to what ha s been observed in many bursts l|Sakamoto et all 
l2007i . 120061 ). Note also that the addition of terms propor- 
tional to Q(r) can not make the spectrum softer because 
they have the same form of energy dependence as the dom- 
inant terms. 

Next we consider the regime where (j'/u)' m > 1. In this 
regime there is no analytical expression for 1F2 functions. 
Therefore we use the asymptotic behaviour of Bessel func- 
tions, K v (x) w y /, ir/2xe~ x , to estimate the integrals in (|B1[) . 
Integration of the left hand side of (|B4|) leads to: 



-(p+i) 



^5/3 (CR; 



—{#(—) 

7m 1 V 2 V ' 



r(-( P + 3),(— )*)_r(-(p + 3),i) 



+ 



(- 



-r 

1 

r(- 



r(f) 

£ 1 I 

a ~ e. 



£ _1_ £ 
4 ~ 6 



11 



5 p 
3' 4 



+ 



6 



E + l. 

4 + 6' 



+ 



for 



> 1 (81) 



The constant coefficient in the last term is ~ 1.354(2 5//3 /(p+ 
2/3) + 2 1//3 /(p+ 10/3)). Using an expansion in series of the 
incomplete F function: 



r(a, x) = e 



x + 



1+7 



(82) 



we see that in this regime at the lowest order the spec- 
trum is a power-law with an index ~ — (1 + p/2) which is 
much steeper (softer) than when uj' /w' m < 1. For energies 
between these two extremes regimes one expect an index 
— (1 +p/2) < a < —2/3. They can be fitted by a power-law, 
but at high u' jui' m a power-law with an exponential cutoff 
gives a slightly better fit due to the presence of the expo- 
nential term in the expansion of the incomplete F functions. 
Nonetheless if the distribution of electrons is truncated at 
high Lorentz factors or is a broken power-law in which the in- 
dex becomes smaller at higher energies, this constraint will 
not be applicable. As ui' /uj' m > 1, the second term in the 
sum and the last term in (|81[1 are smaller than other terms 
considered above. 

The two regimes explained above cover the maxi- 
mum and soft wing of the index range observed by Swift 
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l|Sakamoto et al\ [20071 ') and BATSE but it can not explain 
very hard short bursts with indices larger than the zero or- 
der index —2/3. Moreover, when the observed spectrum is 
extended to high energies, > lMeV, a power-law with an 
exponential cutoff at high energies is a better fit to GRBs 
spectrum. This type of spectra can be obtained if the dis- 
tribution of the Lorentz factor of the accelerated electrons 
has an exponential cutoff at high energies. As we mentioned 
in Sec l2.2l the spectrum of the accelerated electrons in the 
simulations of ultra-relativistic shocks is more sophisticated 
than a simple power-l aw and is closer to a power-law with 
an exponential cutoff ()Spitkov skv 2008), see (|16[1 . Assuming 
for simplicity an exact power-law with exponential cutoff for 
the Lorentz factor distribution of accelerated electrons like 
cq. (120 . the dominant term of the spectrum (|62[) becomes 
proportional to: 



val in which the evolution can be ig nored. Evidently, other 
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with the conditions mentioned in (|20[) . When 7 e <C "/cut the 
exponential term is close to 1 and the behaviour of the spec- 
trum is indistinguishable from a simple power-law. But at 
high energies photons are emitted preferentially by faster 
electrons for which the exponential term is important and 
leads to an exponential cutoff in the spectrum of the syn- 
chrotron emission. To prove such a behaviour, we concen- 
trate on the part of the spectrum for which u' /u' m > 1. In 
this case the integration of (|83[) leads to: 
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For u)' /u' m ^> 1 the second integral on the right hand side 
of (|84[) is small and negligible. The first integral does not 
have analytical solution, but using asymptotic behaviour of 
incomplete T function (|82[l one can conclude that its be- 
haviour at high energies must be close to an exponential 
with negative exponent proportional to uj' . 

In the beginning of this section we mentioned that most 
probably Ar is also somehow energy dependent, specially at 
very low and very high energies. In the latter case it must be 
very small, probably exponentially decreasing with energy. 
This can also contribute to the existence of an exponential 
cutoff at high energies. As for the low energy tail of the 
spectrum the fact that even the residual energy of cooled 
electrons can be enough to emit soft photons means that 
the spectrum in this regime should be much flatter. This is 
consistent with the roughly flat spectrum and shallow tem- 
poral decline in optical an d longer wavelength emissions as 
obser ved in most GRBs (|Roming et al\ 120081 : lOates et al\ 
12008ft . 

In conclusion of this section, we showed that when a 
relativistic shock is modelled in details and realistic distri- 
butions for electrons are considered, the synchrotron emis- 
sion alone can explain different aspects of the time averaged 
spectrum of GRBs as well as their spectrum in a short inter- 



processes such as inverse Compton (IKobavashi et al. 



iPiran et al.ll2008T ) and pair production ( Asano et al 



2007; 



2008) 



contribute to the total emission but most probably are not 
the dominant component at low/intermediate energy bands. 
Nonetheless their contrib ution should be mo re important at 
GeV and higher energies (I Asano et aLll2008T ). 



4 SUMMARY 

We presented a formulation of the relativistic shocks and 
synchrotron emission that includes more details than the 
dominant term considered in the previous calculations. Al- 
though we consider spherical shells, most of our results are 
valid also for non-spherical collimated jets as long as the col- 
limation angle due to the relativistic boost is smaller than 
collimation angle. 

We showed that the lags between light curves at differ- 
ent energies exist in the dominant order and are not only 
due to the high latitude emissions which are negligible for 
ultra-relativistic ejecta. The main reason for such behaviour 
is the evolution of electric and magnetic fields as well as the 
evolution of the emitting region which can be in addition en- 
ergy dependent. This fact is more evident in the simulations 
presented in Paper II. Despite the absence of high latitude 
terms in our simulations, the presence of lags between the 
light curves of different energy bands is evident. For the same 
reasons the change in the slope of the light curves - what 
is called breaks - are also energy dependent. This explains 
chromatic breaks of the GRBs detected by Swift . 

The two phenomenological models we considered for the 
evolution of the active region are physically motivated, but 
do not have rigorous support from microphysics of the shock. 
Nonetheless, they can be easily replaced if future simula- 
tions of Fermi processes lead to a better estimation of the 
size of the region in which electric and magnetic fields are 
formed and particles are accelerated and dissipated. Pres- 
ence of an external magnetic field in the environment of the 
candidates for central engine of GRBs is very plausible. The 
formulation present here does not include such a possibility, 
but an external magnetic field can be added to (|23|l . The 
modification of the evolution equation of j3' and the flux is 
straightforward . 

Many other details such as the effect of metalicity of 
the ejecta and surrounding material both on the low energy 
emission and absorption is not considered in this work. We 
have also neglected synchrotron self-abortion. It only affects 
the low energy bands, nonetheless in hard bursts even optical 
emission can be affected by self absorption. We leave the 
study of this issue, the effect of ionization on the emissions, 
and the thermalization of shocked material to future works. 



APPENDIX A: ANALYTICAL 
APPROXIMATION OF (5' (1) 

The integral in (I3T[) can not be taken analytically. An ana- 
lytical approximation is however useful for investigating the 
effect and importance of the radiation in the dynamics of 
the shell collision. Here we consider a few approximations 
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for various range of T> and 77 that determine respectively the 
strength of the shock and the radiation. 

If 77 = T = 0, the integral in (|37[) can be calculated 
analytically: 
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When 77/ (3— ft) is large, as the value of y in the integrand of 
(|37[) is always less than 1 we can formally expand the term 
(1 — y 2 ) 3 ~ T ^ 2 and integrate term by term: 
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Here 2*1 is the hypergeometric function and can be ex- 
panded as a polynomial of its last argument. When the latter 
is less than 1 and the power of the terms in the polynomial 
are positive, they converge rapidly to zero. When the last 
argument in 2 F\ is larger than one an analytical extension 
of this functio n with negative power in the polynomial ex- 
pansion exists (|Gradshtevn fc Rvzh ik 1980). Therefore, the 
dominant term is always /3[ j in front of each 2 F\ term. For 
small 77 this approximation is not valid. In this case we ex- 
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In the same way we find the following expression for 
Mi ){r') for the quasi steady model and large 77/(3 — ft): 
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For the case of a small 77/(3 — ft) one can use (|A4[I and make 
an expression analogue to (|A5|I . thus we do not repeat the 
details here. We note that when t — and S — > 00, i.e. for a 
constant Ar', (IA3|) and (|A5[) are equal as expected. When 
5^0, the formation of the active region is very slow. In this 
case the energy loss by radiation becomes negligible and j3 
decreases only due to the increasing accumulated mass. Note 
also that A4(o)(t"o) = as expected. If the active region 
varies according to ([29}, the expression for M(p)(r') includes 
only the term 77 — > T] + 5 in (|A5I) with a positive sign in front. 

Although equations (|A3|) . (|A5|) . and (|A4|) seem quite so- 
phisticated, due to the polynomial representation of 2 Fi, 
only the few dominant terms are of real interest to us. In 
most cases we are only interested in the dominant power-law 
component. However, having expressions beyond the domi- 
nant power permits to go much further and calculate quan- 
tities such as lags that in a simple power-law approximation 
can not be determined. 

When the kinetic energy of the fast shell does not 
change significantly ft )( r ') ~ ft and M(p)(r) — > 0. This 
case happens when the radiation has a negligible effect on 
the kinematic of the shock. Assuming a strong shock i.e. 
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P' < 3D, we can use the definition of %F\ to investigate the 
behaviour of Mro)(. r ) at lowest order. For the dynamically 
driven active region with large r\ and a relatively soft shock, 
when we expand 2-F1 terms in (|A3|I up to first order, M (0) (r) 
becomes: 



APPENDIX B: SPECTRUM AND LAGS FOR 
POWER-LAW AND EXPONENTIAL 
ELECTRON DISTRIBUTIONS 

Using (|13[) to (|15[) . after integration over 7 e we obtain the 
following expression for the spectrum: 
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Similarly, for small r\ we use (|A4[I and obtain: 
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For a quasi static active region and the same shock condi- 
tions as in (|A5f) the first order _M(o)(f) is: 
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where cj^ = Wcc7m i s the minimum characteristic frequency 
of electrons. The double integrals in the first and last terms 
of (|B1|) do not have a simple analytical expression. An ap- 
proximation can be obtained using an integral form of K v 
Bessel function: 
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The last integral is small and can be neglected. Therefore: 
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If the electron distribution is exponential i.e. n' e ('fe) = 
N e exp(—ccy e ), terms with 1F2 must be replace by terms pro- 
portional to Mejier's G-functions. These functions do not 
have simple analytical presentations. Therefore only asymp- 
totic behaviour of power spectrum is described in Sec. 12.41 
In the same way we can determine an analytical expres- 
sion for lags when accelerated electrons have a power-law 
distribution using (J65J) to J75]) and (|BT|) to (fB5l) . In fact, it 
is easy to see that the functions P(ro,w'o) and Q(ro,u}' ) are 
both functionals of Ti which contains 7 e dependent terms 
similar to the spectrum dP/coduj. The function Tt includes 
Bessel and hypergeometric functions and their derivatives. 
Therefore we can use (|B1[) to (|B5I) to determine them. As the 
derivatives of hypergeometric functions \F2{a, /3, 7; z) are 
also hypergeometric, lags can be expressed as a sum of 1F2 
functions. The calculation is straightforward but long and 
laborious and we do not present details here because their 
complexity does not permit to investigate their properties 
and numerical calculation is needed. 
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